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Abstract 

We study the perturbative QCD description of the HERA measurements of F2(x, Q 2 ) 
using a gluon distribution that is obtained from an evolution incorporating angular or- 
dering of the gluon emissions, and which embodies both GLAP and BFKL dynamics. 
We compare the predictions with recent HERA data for F2. We present estimates of the 
charm component F£(x, Q 2 ) and of Fl(x, Q 2 )- 



1 On leave from Henryk Niewodniczariski Institute of Nuclear Physics, 31-342 Krakow, Poland. 



Deep inelastic electron-proton scattering experiments at HERA have measured the structure 
function F 2 (x,Q 2 ) in the previously unexplored small x regime, x < 10~ 3 . The values of F 2 
are found to rise rapidly with decreasing x . These measurements have stimulated much 
theoretical activity and the small x behaviour of F 2 has been interpreted using perturbative 
QCD from several different viewpoints. The interpretation is complicated by the need to provide 
non-perturbative input. 

In fact the present data for F 2 can be well described by traditional Altarelli-Parisi (or GLAP) 
evolution in the next-to-leading order approximation. The data imply a steep gluon (that is a 
gluon density which increases as x decreases) even at low Q 2 values. We may input the steep 
x behaviour directly into the starting distributions at some input scale, say Ql = 4 GeV 2 0, 
or alternatively we may generate it from "non-singular" or "flat" x distributions at some low 
scale, such as Ql = 0.3 GeV 2 @] or Ql = 1 GeV 2 ||, chosen so that the evolution length 

rQ 2 da 2 
£(Q 2 ,Q 2 )= / « s ( 9 2 )^ 
JQ q 

is sufficiently long, where as = 3as/7r. GLAP evolution amounts to the resummation of the 
leading (and next-to-leading) logQ 2 terms. At small x and large Q 2 IQ\ it generates a steep 
double leading logarithmic (DLL) behaviour] 2 ] of the form exp(2[£(Qo, Q 2 ) log(l/x)] 2 ). Despite 
the apparent success of the GLAP approach, it is not the only way of generating a steep gluon 
compatible with the HERA data. 

At sufficiently small x we must also resum the as logl/x terms, unaccompanied by logQ 2 . 
This is accomplished by the BFKL equation. It generates a singular x L behaviour of the 
unintegrated gluon distribution, f(x,k^), where Xl = ag41n2 for fixed as- If a reasonable 
assumption is made to introduce the running of 0:5 then the numerical solution of the BFKL 
equation again yields an x~ x behaviour but with A ~ 0.5 0. Using the /^-factorization theorem 
the behaviour f(x, k^) ~ x~ x feeds through into F 2 (and into Ft)- To be precise we have 

F t (x, Q 2 ) = J^f £ % fW> 4) ^ b ° x 4, Q 2 ) + if (1) 

with i — 2,L, where Ff ~ Fi(x,Q 2 ) at large x, but is a slowly varying function of x and Q 2 
at small x. The convolution in (p]) is diagrammatically displayed in Fig. 1. _p. box includes both 
the quark box and crossed box contributions which originate from virtual photon- virtual gluon 
qq production. For the cc component we take the quark mass to be m c = 1.4 (or 1.7) GeV. At 
sufficiently small x the x~ x BFKL behaviour overrides the DLL form. To find precisely where 
this will happen requires a unified BFKL/GLAP formalism, as well as knowledge of the yet 
unknown next-to-leading logl/x contributions. 

2 This form increases with decreasing x faster than any power of logl/a; but slower than any power of 1/x. 
A choice of singular starting distributions, xg, xq sc & ~ x~ x with A > 0, would therefore eventually override the 
DLL behaviour. 
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One approach which unifies the GLAP and BFKL formalisms is to recast the leading twist 
part of the BFKL fey-factorization formula into collinear form in which the splitting and co- 
efficient functions acquire the higher order log 1/x contributions J7|. Several interesting phe- 
nomenological studies have developed from this formalism R IS], ITUf . Actually the procedure is 



to use the BFKL equation for fixed as to obtain as/co power series expansions of the anomalous 
dimensions and coefficient functions, and then to let as run. Here uj is the moment index. In 
this way we retain the simplicity of GLAP evolution, but with splitting and coefficient functions 
which incorporate BFKL resummationsf]. 

We can see that the basic quantity is the unintegrated gluon distribution /, which corre- 
sponds to the sum of (effective) gluon ladder diagrams. The unintegrated gluon satisfies the 
BFKL equation 



f(x,k 2 T ) 



1 dz f dk% 



Jx Z J 



U2 



f(z,k'*)-f(z,k 2 T ) 

I h' 2 - h 2 I 



+ 



f{z,k\ 



(4A# + A£)3_ 



(2) 



where a~s = 3as/7r. We notice the potential collinear singularity at k'^ = 0. However, provided 
the driving term f° is chosen to vanish at k\ = 0, the structure of the equation guarantees 
that it is free of collinear singularities. This is a natural wayQ to regulate the singularity. It 
can be linked directly with the Qo scheme advocated by Ciafaloni ITD]. Since we stay in four 



dimensions we avoid factors which are characteristic of minimal subtraction (MS) schemes. 

Here we work with the fey-factorization formula, (|l]), and do not reduce the equation to 
collinear form. This allows us to study the effect of replacing the BFKL gluon with the gluon 
obtained from the CCFM equation a unified equation which embodies both the BFKL 
equation at small x and GLAP evolution at large x. The CCFM equation is based on the 
coherent radiation of gluons, which leads to an angular ordering of the emitted gluons. 

Since the next-to- leading log 1/x contributions are not yet known, the introduction of run- 
ning as into the BFKL equation]] is, of necessity, subject to assumption. The most reasonable 
procedure is to take asik?) in (||), so that the equation is compatible with the double- leading- 
logarithm limit of GLAP evolution. This prescription, however, generates a solution which 
differs from that obtained when the (fixed as) BFKL equation is first reduced to collinear form 
and then as is allowed to run JHJ. In other words the introduction of asik 2 -) in (fj) gives a 
different solution to that obtained by allowing as to run in the leading-twist collinear solution 
of the BFKL equation. Formally, however, the difference between the two approaches can be 
attributed to non-leading In 1/x effects. For instance in the region 



«s(Q 2 



< 



ois{Q. 



< 



CO 



CO 



41n2 



3 There is freedom in assigning the BFKL resummations to the coefficient functions, splitting functions and 
the starting distributions. Various factorization schemes have been proposed [O, |l2|, |l3| . 

In fact the vanishing of the inhomogeneous term is ensured by the colour neutrality of the probed hadron. 
5 We thank M. Ciafaloni and L. Lipatov for valuable discussions on this point. 
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both methods give the same result up to these non-leading effects. However, the singularity 
structure in the moment or uj plane is different in the two formulations. The collinear reduction 
of the BFKL equation, leading to conventional evolution from, say, Qq to Q 2 , contains a branch 
point singularity at uj = ul(Qq) = a s(Qq) 4 In 2 or, if this leading singularity is absorbed in the 
starting distributions, at uj = ojl{Q 2 )- On the other hand the solution with running a$ directly 
incorporated into the BFKL equation does not contain the branch point singularity but rather 
it has (an infinite number of) poles in the uj plane. The leading pole is well separated from the 
others with a position uj v pa aasibk^) where a and b are constants (b pa 7) and k delimits the 
infrared region. It turns out that uj p < uj^k^) |R| . 



Following ref. ||16|| , we now implement angular ordering of the gluon emissions. The uninte- 



grated gluon distribution is then a solution of the CCFM equation^] JT6| rather than the BFKL 



equation. In the small x region the CCFM equation may be approximated by 

f(x,k 2 T ,Q 2 ) = f(x,k 2 T ,Q 2 )+a s k 2 T £^ J ^Q(Q- zq )A R (z,q,k T )±f(~,k' 2 ,q^ (3) 

with k' T = krp + q, see eq. (18) of ref. JT7|. A R represents the virtual corrections which screen the 
l/z singularity and the theta function imposes the angular ordering on the real emissions. Eq. 
(11) of ref. |17| gives the explicit expression for Ar. We note that the solution / depends on an 



additional scale Q that is required to specify the maximum angle of gluon emission (which turns 
out to be essentially the scale k of the probe, see Fig. 1). This equation has recently been solved 
numerically and the resulting gluon distribution has been compared with that obtained from 
the BFKL equation [17| . As anticipated, the angular ordering constraint suppresses the CCFM 
gluon at the lower Q = k 2 values. If we replace the angular-ordering constraint Q(Q — zq) 
by Q(Q — q) and set Ar = 1 then we obtain an equation which becomes equivalent to the 
Altarelli-Parisi (GLAP) equation in the double- leading- logarithm approximation (DLL A). 

The small x approximation of the CCFM equation that we have used (see WW l amounts 
to setting the Sudakov form factor A$ = 1 and to approximating the gluon-gluon splitting 
function by its singular term as z — > 0, that is P gg ~ 6/z. As represents the virtual corrections 
which cancel the singularities at z — 1. We account for the remaining finite terms in P gg by 
multiplying the solution f(x, k\, Q ) by the factor 



exp 

where A is defined by 



A f Vstf) f) (4) 



r 1 6 

^ z" P gg {z)dz ~ - - 6A (5) 

That is A = (33 + 2?2y)/36, where the number of active flavours nj = 4. 

Fig. 2 compares the CCFM and DLLA predictions for F 2 with the recent HERA mea- 
surements 0, H- The predictions are obtained by first determining the gluon distribution 

6 The CCFM equation incorporates part of the non-leading logl/a; contributions. 
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f(x, fcfi, Q 2 ) by iteration of (|3|)[] in the domain k\ > k\ = 1 GeV 2 starting from a "flat" driving 
term of the form 3(1 — x) 5 exp(— k^/k 2 ) with k 2 = 1 GeV 2 , that is exactly as in ref. [ I7| - We 



correct for the small x approximation by multiplying the gluon distribution / by the factor 
shown in (Q) and then predict F 2 from the /^-factorization formula ([!]) with an infrared cut-off, 
k\ > k%. For F 2 5 we use the value of F 2 (x,Q 2 ) obtained from the MRS(A') set of partons || 
at x = 0.1, and extrapolate below 0.1 assuming the normal x -0 - 08 "soft" behaviour. 

From Fig. 2 we see that the CCFM and DLLA predictions coincide at large x, as indeed 
they should. The two schemes start to differ at small x and Fig. 2 indicates the value of x 
at which the resummation effects become important. It is evident that once a background is 
added to the small x behaviour predicted by the CCFM equation then a good description of the 
HERA data is obtained. The prediction lies between the GRV and MRS (A') values. It should 
be noted that the CCFM calculation is not a fit to the HERA data, but simply a solution of 
the evolution equation incorporating angular ordering. The rise of the gluon, and hence of F 2 , 
is generated by the evolution equation and hence is within the domain of perturbative QCD. 
Of course the perturbative QCD prediction is not absolute. The normalisation depends on the 
choice of A; 2 , which delimits the infrared region, and also on the choice of the driving term. 
Also the normalisation depends on the choice of the lower limit of integration in (|4]). Here 
we take this to be Q = 1 GeV 2 . Recall that the correction factor (£|), and hence Q , only 
occurs because we solve a simplified form of the CCFM equation appropriate to the small x 
region. In summary there is some freedom in the normalisation of F 2 , though the prediction of 
the shape of the x dependence is characteristic of the CCFM equation. It is encouraging that 
the physically reasonable choice fc 2 = Q = 1 GeV 2 gives such a satisfactory description of the 
HERA data. 

For completeness we show in Figs. 3 and 4 respectively the predictions for the longitudinal 
structure function Fl(x, Q 2 ) and for the charm component of F 2 , which we denote by F£(x, Q 2 ). 
In each case we specify the background or "soft" contribution Ff at x — 0.1 to be given by 
the MRS (A') predictions, and extrapolate below 0.1 using the x -0 ' 08 "soft" behaviour. For the 
charm component F% the CCFM predictions for m c = 1.4 and 1.7 GeV are shown; the argument 
of the running coupling is taken to be k 2 + m 2 (where k 2 is shown in Fig. 1). The predictions 
of Fl and F^ obtained from GRV and MRS (A') partons are also shown in Figs. 3 and 4. 

Since the CCFM values of F 2 agree with the HERA data, we can regard the charm com- 
ponent F£ as an absolute prediction. The charm component of the MRS (A') partons has been 



fixed to be in agreement with the EMC measurements [ITS]] of F 2 C which lie in the region x ~ 0.1. 
Indeed we see these data barely extend into the kinematic region shown in Fig. 4. It will be 
particularly informative to have measurements of F£ at HERA in the small x regime where 
resummation effects are expected to occur. 

In summary we have shown that it is possible to obtain a good description of the HERA 
measurements of F 2 from the solution of a unified evolution equation based on the angular 



7 The DLLA prediction is obtained as defined above, that is by taking A/? = 1, and z = 1 in the O function 
§)■ 



4 



ordering of the emitted gluons. The gluon distribution f(x,k?,Q ) was obtained by iteration 
starting from a driving term of the form 3(1 — x) 5 exp^—k^/k^), and the structure function F 2 
was then determined via the /^-factorization formula F 2 = f <8> -F 2 box . The steepness of the 
gluon, and of F 2 , with decreasing x, is generated by the evolution equation. In this way we 
identified the regime where the ln(l/x) resummations become important. 

However, our treatment is only a first step. There are several reasons why it may over- 
estimate the rise, particularly at low Q 2 . First we have to find a realistic way to impose 
energy-momentum conservation of the emitted gluons. Second, we have ignored gluon shad- 
owing corrections. These are expected to be small in the HERA regime, as evidenced by the 
persistent rise of the F 2 data with decreasing x for Q 2 as low as Q 2 = 2 GeV 2 . Last, but not 
least, the full next-to- leading ln(l/x) contribution is unknown at present. This is needed to 
check the prescription for the running of as and to specify the scale dependence. 

Clearly the agreement of our CCFM predictions with the small x measurements of F 2 do not 
imply angular ordering effects have been firmly established. GLAP and BFKL evolution can 
give an equally good description. There are two characteristic features of the gluon distribution 
f(x,k^,Q 2 ) obtained from an evolution equation which includes a resummation of \n(l/x) 
terms. Namely a steep rise of / with decreasing x which is accompanied by a diffusion in 
ln/cf.. F 2 measures only the rise. A distinctive test will involve both features. For this we 
need to explore final state processes such as deep inelastic events containing an identified 
energetic forward jet. Here we have focused on F 2 and obtained predictions based on angular- 
ordered evolution which embodies both BFKL and GLAP resummations. Moreover, we have 
also presented values for the charm component F 2 and the longitudinal structure function Fl. 
Valuable insight into the properties of the gluon can be obtained from the measurement of 
these structure functions at HERA. 
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Figure Captions 

Fig. 1 Pictorial representation of the fc^-factorization formula, that is of the convolution F t = 
/<S>.Fj box contained in eq. (]l|) with % = 2, L. f(x', kj) is the unintegrated gluon distribution 
and Fj box is the off-shell gluon structure function, which at lowest order is determined by 
the quark box (and "crossed" box) contributions. The integration variables, x' and k^, 
are respectively the longitudinal fraction of the proton's momentum and the transverse 
momentum carried by the gluon which dissociates into the qq pair. 

Fig. 2 A comparison of the HERA measurements of F 2 [|1|, |2| with the predictions obtained 
from the /^-factorization formula (|I|) using for the unintegrated gluon distribution / the 
solutions of the CCFM equation (continuous curve), and the DLL-approximation (dot- 
dashed curve) of this equation. We also show the values of F 2 obtained from collinear 
factorization using the MRS (A') || and GRV [|j] partons. 

Fig. 3 The continuous curve is the prediction for the longitudinal structure function of the 
proton, Fl(x,Q 2 ) obtained by solving the CCFM equation for the gluon and using the 
/cT-factorization formula ([I]). The values of Fl obtained from GRV Q and MRS (A') || 
partons are also shown. The GRV prediction includes a charm component only at leading 
order. 

Fig. 4 The continuous curves are the predictions for the charm component F% of the proton 
structure function F 2 obtained by solving the CCFM equation for the gluon and using 
the /cT-factorization formula with m c = 1.4 (upper) and 1.7 GeV (lower curve). The 
values of F 2 obtained from GRV at leading order f|] and from MRS (A') || partons are 
also shown. The next-to-leading order GRV prediction lies below the leading order result; 
for example at Q 2 = 15 GeV 2 and x = 10~ 4 it is shown by a small cross. Also shown are 
EMC data |L8| at adjacent (Q 2 ) values, assuming that the c —>■ fi + X branching ratio is 
8%. 
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Fig. 2 
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